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Understanding the effect of inhomogeneity on the charge regulation and dielectric 

m ; 

properties, and how it depends on the conformational characteristics of the macro- 

"xT ■ 

\& , molecules is a long-standing problem. In order to address this problem, we have 

o 

developed a field-theory to study charge regulation and local dielectric function in 
planar polyelectrolyte brushes. The theory is used to study a polyacid brush, which 



is comprised of chains end-grafted at the solid-fluid interface, in equilibrium with a 
bulk solution containing monovalent salt ions, solvent molecules and pH controlling 
acid. In particular, we focus on the effects of the concentration of added salt and pH 
of the bulk in determining the local charge and dielectric function. Our theoretical 
investigations reveal that the dipole moment of the ion-pairs formed as a result of 
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counterion adsorption on the chain backbones play a key role in affecting the local 
dielectric function. For polyelectrolytes made of monomers having dipole moments 
lower than the solvent molecules, dielectric decrement is predicted inside the brush 
region. However, the formation of ion-pairs (due to adsorption of counterions coming 
from the dissociation of added salt) more polar than the solvent molecules is shown 
to increase the magnitude of the dielectric function with respect to its bulk value. 
Furthermore, an increase in the bulk salt concentration is shown to increase the local 
charge inside the brush region. 



INTRODUCTION 



Polyelectrolyte brushes[l-6] have been used extensively in the development of stimuli- 
sensitive materials due to their ultra-sensitive response to different stimuli such as pH, tem- 
perature, solvent etc. A fundamental understanding of the effect of each of these stimuli on 
the brush properties is highly desirable for a number of technologies. Extensive research jj- 



22J aimed at elucidating the effect of different experimental variables on the response of 



arushes has been conducted. However, in general, dielectric properties of 



23 



271 ] elude a clear understanding and are rarely studied using rigorous 



polyelectrolyte 
polyelectrolytes 
theoretical tools. 

Classic literature on the dielectric function 



28 



33| mainly deals with the small molecules 



randomly distributed in space, constituting a spatially homogeneous medium. For a spatially 
uniform medium, in the absence of explicit charges, the dielectric function also turns out 
to be invariant in space and hence, the words "dielectric constant" are used in place of 

- l36| that 
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"dielectric function" . For a single charge in a polar medium, it has been shown 
the dielectric function depends on the distance from the charge. The dielectric function 
is lowest next to the charge and monotonically approaches the dielectric constant of the 
medium. This effect is known as the dielectric saturation effect 28j resulting from the strong 
electric field next to the charge, which vanishes with an increase in the distance from the 
charge, reaching zero in the bulk. 



The dielectric saturation effect has been shown [28j to be responsible for the dielectric 
decrement when trace amounts of electrolytes like sodium chloride are added to polar sol- 
vents such as water. In particular, a linear dependence of the dielectric decrement on the 
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salt concentration was proposed by Debye 



The linear relation between the dielectric 



decrement and the salt concentration was also observed in experiments [28]. However, the 
situation at large concentrations of salt is far more complex due to the interplay between 
association/ion-pairing and short-range interactions. 

Like electrolytes, polyelectrolytic systems are inherently inhomogeneous. Understanding 
the effect of inhomogeneity on the dielectric properties of the macromolecules is a long- 
standing problem. To the best of our knowledge, there is no comprehensive theory for the 



local dielectric function of polyelectrolytes, despite a large body of experimental work[23N26|. 



which has shown that addition o 
salt leads to dielectric increment 



: polyelectrolytes to a solution containing polar solvent and 



23j-|26(. Minakata et al. |27| attributed the overall dielectric 



increment to the fluctuations of the bound counterion charge at different sites. However, 



the theory by Minakata et al. 



271 ] ignores the conformational degrees of freedom and deals 



with rod-like polyelectrolytes in an external electric field. 

From the above discussion, it is clear that the dielectric function depends on the number 
of charges and their distribution. This dependence of the dielectric function adds another 
complexity in the case of polyelectrolytes. The complexity arises because of the tendency of 



polyelectrolytes to regulate their charges 



37H4l| on the basis of the local environment such 



as local pH, salt cocentration, etc. In other words, one has to compute the local charge 
and the dielectric function in a self-consistent manner while minimizing the free energy. In 
this work, we have developed a quantitative description of the charge regulation and local 
dielectric function for polyelectrolytes. The theoretical formalism is quite general and we 
use it to study a planar polyelectrolyte brush in equilibrium with a solution bath containing 
monovalent salt ions, solvent molecules and an acid to control the pH of the bath. This 
particular set-up serves as a model system and is relevant for a number of experimental 
studies. Typically, solvent molecules are treated as a continumm without having any struc- 
ture. Here, we treat each solvent molecule as having a permanent dipole moment of its 
own, which is a reasonable model representation of a number of polar molecules. For a 
quantitative analysis of the dielectric function, which must depend on the local properties 
such as density, electrostatic potential and electric field, we have developed a field-theory 
taking into account the dipolar interactions between charged and uncharged components. 

This paper is organized as follows: the general formalism is presented in section [Til and 
details about the application of the formalism to the polyelectrolyte brushes are presented 
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in section II III Numerical results are presented in section IIVI and section |V] contains our 
conclusions. 



II. THEORY 

We consider a planar polyelectrolyte brush formed by n mono-disperse flexible polyacid 
chains (such as poly(methacrylic acid)), each having N Kuhn segments of length I. The 
chains are assumed to be uniformly grafted onto an uncharged substrate so that the grafting 
density (defined as the number of chains per square nanometer) is a (see Fig. [1]). For 



the field theoretical analysis 



42 



43| described in this work, each polyelectrolyte chain is 
represented by a continuous curve of length Nl, and an arc variable t is used to represent 
any segment along the backbone so that t £ [0, N]. t = corresponds to the grafted end 
and t = N represents the other end. To keep track of different grafted chains, subscript 
a is used so that t a represents the contour variable along the backbone of a th chain. In 
the following, we use the notation H a (t a ) to represent the position vector for a particular 
segment, t a , along the a th chain. 

Counterion adsorption on the chains is studied using a two-state model. In this model, the 
entire population of counterions is divided based on whether they are "free" or "adsorbed" . 
As the label implies, the "free" counterions are free to sample the whole space, while the 
"adsorbed" counterions have less translational degrees of freedom and are bound to the 
chains. In other words, each segment of a chain can be in either uncharged or charged 
state. The adsorbed counterions/uncharged segments are treated as "point" dipoles. For 
each uncharged segment t a along the a th chain, an electric dipole of moment (in units of 
electronic charge, e) p a (t Q ) is assigned. Similarly, each solvent molecule is assigned an 
electric dipole moment and we use the notation to represent the dipole moment of the 
k th solvent molecule. Interactions between the charged and uncharged species are taken into 
account by ion-dipole potentials. 

In this work, we focus on the effects of "permanent dipole moments" and don't allow 
variations in the magnitude of the electric dipole moments. However, we allow the dipoles 
to be aligned depending on the local electric field. The physical origin of these electric dipoles 
is the difference in electronegativities of the atoms constituting the dissociable group along 
the backbone. For example, acetic acid has a dipole moment of 1.7 Debye[44]. It is well- 
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known 28M31| that finite polarizability of molecules leads to an enhancemnet of net electric 
dipole moments. We plan to include these effects in future work. 




t = 





t = N 



FIG. 1: Schematic of a planar polyelectrolyte brush containing monodisperse chains in equilibrium 
with an electrolyte solution (not shown explicitly here) . The polyelectrolyte chains are modeled as 
continuous curves containing charges (represented by red dots) and electric dipoles (represented by 
arrows) resulting from counterion adsorption. The pH regulating electrolyte solution is assumed 
to contain two kinds of counterions resulting into two kinds of dipoles on the chains. 



We assume that the brush is in equilibrium with a bulk solution whose pH is controlled 
by the use of a buffer containing an acid and the precise concentrations of protons, H + , 
cations from the added salt, B + , and anions, A~, which result from the dissociation of the 
acid, water molecules (protons and hydroxyl ions) and the salt, are known. For the sake 
of simplicity, we assume that the anions coming from the acid, water molecules and the 
salt are indistinguishable. Generalization of the theory presented here to take into account 
the specificity of anions is straightforward. However, we consider the specificity of cations 
(counterions for the brush) to study the effects of dipole moments of different kinds of ion- 
pairs and competitive counterion adsorption. Physically, H + originates from the dissociation 
of acidic sites on the polyelectrolytes, protonation of water molecules and the acid (for pH 
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< 7). In this paper, we use the word "bulk" quite frequently while referring to the region 
in space, where electrostatic potential and densities are spatially independent. One should 
not confuse the region just outside the brush region as the "bulk", because the electrostatic 
potential may be inhomogeneous even in monomer-free regions. 

The theory is developed in the Canonical ensemble and the equilibrium with the bulk so- 
lution is considered by equating the chemical potentials of different components in the inho- 
mogeneous brush and the homogeneous bulk solution. For the theoretical treatment carried 
out in the Canonical ensemble, we assume that there are n 1 ions of kind 7 = H + , B + , A~ . 
The number of counterions in the "free" and "adsorbed" states are represented by super- 
scripts / and a, respectively. Overall, the system is electroneutral and local electroneutrality 
is assumed for the bulk solution. To keep track of different kinds of ions, we define Z 7 as 
the valency (with sign) of the charged species 7 = H + , B + , A~ ,p, where the subscript p 
represents the polymers. For monovalent ions, \Z^\ = 1 and for a polyacid, Z p = — 1. 

Following the mathematical procedure presented in Appendices A and B, we can write 
the partition function for the polyelectrolyte brushes as 



Z = \j D \p p ] J D [w p ] J D [w s ] J D exp 



H 



k R T knT 



(1) 



where 



k B T 



"T^F + o l w ppPponN + w ss p so n s 



In 
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(2) 



and F a contains contributions coming from the adsorbed counterions. Explicitly, it is given 
by (see Appendix B) 



k B T 



n a B+ lnK B + - (nN - n a H +) \nK H + - In 



n a H+ \n a B+ \{nN — n a H+ — n 



+ ln 



n f H+ \n f B+ \n s \ 



- (nN + n s )ln47r 



(3) 



In writing Eq. [3j we have defined an equilibrium constant, K H +, for the dissociation 



45| K 



H+ 



of acidic groups on the polyelectrolyte chains by the thermodynamic relation 
exp[— {n° coo - + P° H + — Pcooh}} where \i is the chemical potential (in units of k B T, the 
Boltzmann constant times the absolute temperature) and superscript o represents the limit of 
infinite dilution. Similarly, we have defined another equilibrium constant K B + for the binding 
of the counterions from the salt by the relation K B + = exp[— {f^coo- + P°b+ ~ ^000-3+}]- 
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Also, 
H 



k B T 



\,,, I drp p (v)p 
1 



1 - 



Pp( r ) 



Ppo 



+ i / dr [V> p (r) - w p (r)] p p (r) 



8ttZ 



Bo 



dr^(r)V^(r) 



y drp so 


J Pp( r )~ 




u> s (r) + In 


Ppo _ 







sin (p 8 |V r ^(r)|) 
p s |V r ^(r)| 



In Qy {ip} - n A - In Q A - {^} - n s In Q s {w s } - ^ Q pa {w p } (4) 

y=H+,B+ 



a=l 



Here, Q 7 , Q<j and Q pQ are the normalized partition functions for a small ion of type 7, solvent 
molecule, and a th grafted chain, respectively, given by 



Q 7 M = drexp[-iZ 7 ^(r)] for 7 = H + , B + , A~ 
Q s {w s } = drexp [-iw s (r)) 



(5) 
(6) 



and 

Qpa {^p} 



D [R Q (t a )]exp 



3 

2P 



N 



dt r 



dH, a (t a t 



'R a (0)=r Q 

where r a is the position vector for the grafted end of the a th chain. 



JV 

— 1 I dt a w p (R, a ) 
'0 



(7) 



In these equations, ip(r) is the collective field introduced to decouple electrostatic in- 
teractions and it is the equivalent to the electrostatic potential. w p (r),w s (r) are the fields 
introduced to decouple short range interactions modeled by the Edwards' delta functional 



approach 



41 



431 ] . p p (r) is the collective density variable. The Edwards' approach for using a 
delta function potential to model short range interactions along with the local incompress- 
ibility constraint naturally leads to appearance of the parameter Xps in the theory, defined 
in terms of excluded volume parameters w pp , w ss and w ps in Appendix B. 

Also, I Bo = e 2 /e kBT is the Bjerrum length in vacuum, where e is the charge on an 
electron and e D being the permittivity of vacuum. Q is the total volume of the system and 
p s = |pfc| is the magnitude of the dipole moment of each solvent molecule. p po and p so are the 
densities of pure polymers and solvent, respectively, and S is a normalizing factor defined 
in the Appendix B. The function ip p appearing in Eq. H]is given by (see Appendix B) 



exp [-# P (r)] 



[1 - H + - Pb+) exp [-iZ p ip(r)) + /3 y 

-y'=H+,B+ 



sin (py 




py 


V P V(r)| 



(8) 
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where /3y , 7' = H + , B + are the variational parameters characterizing the extent of counterion 
adsorption of type 7' on the polyelectrolyte chains. These parameters need to be determined 
by the minimization of the free energy of the system. Furthermore, py is the magnitude of 
the dipole moment of an ion-pair formed by the adsorption of counterion of type 7'. 

In the next section, we present the implications of the dipolar interactions on the proper- 
ties of planar polyelectrolyte brushes in equilibrium with a solution containing monovalent 
salt. 



III. PLANAR POLYELECTROLYTE BRUSH IN EQUILIBRIUM WITH A 

SOLUTION BATH 



A. Saddle-point approximation 



An inhomogenous polyelectrolyte brush in equilibrium with a salty solution can be stud- 
ied using the saddle-point approximation 43] . The approximation evaluates the functional 
integrals over the fields by the value of the integrand at the saddle-point. Optimizing the 
Hamiltonian given by Eq. H]with respect to w p ,w s ,p p and if), respectively, we obtain 

5\nQ pa {w*} 



Ppo 



a=l 

n s exp 



<H( r ) 
-<(r)] 



/ dr exp [— w*(r)] 



(9) 
(10) 



W p [T) 



Pso 
Ppo 



w*(r) + In 



sinh (p s |V r ^>*(r)|) 
P,|V r ^(r)| 



+ XpsP s 



2p*(r) 

Ppo 



V r ^*(r) 



-An 



z YP ;,(v) + z A - P \.{v) + z p /3*(t) p ;(v) 

■y'=H+,B+ 



(12) 



w*{r),iw s [r t 



w*(T),iip(r) = ip*(r),iip p (r) 



43] at the saddle point. Also, the 



Here, we have used the notation iw p (r) 
ipp(r) because the collective fields are purely imaginary 
collective density variables such as p p (r) at the saddle point are as p*(r). p*(r) in Eq. 
is the number density of small ions of type 7 = H + , B + , A", given by 

n 7 exp [— Zy0*(r)] 



P 7 l r J 



J dr exp [— Z-yip*(r)} 



(13) 
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where n 7 = n , for 7 = H + , B + . In other words, only the "free" counterions appear in Eq. 
IT2l Furthermore, the effective Bjerrum length in Eq. [12] is given by 



h{p* p ^*} 



-i- + 4tt J J2 PyPy(r)p;(r)L by |V r ^(r) 
lflo (y=h+,b+ 



+PsPso 



1 - 



P|W 

Ppo 



L[p s |VW>*(r) 



(14) 



where = L(x)/x so that L(x) = cothx — 1/x is the Langevin function. Furthermore 

(l-(3 H+ -/3 B+ )exp {-Z p r(r)] 



^(r) 



exp [-^* 



sinh(py I V r ^*(r) |) 



exp [-^*(r)] Py|V r ^*(r) 



(15) 
(16) 



so that if)* is given by Eq. [S]at the saddle point with the notation iip p (r) = Explicitly, 
it is given by 

exp [--0*( 



;i - Ph+ - Pb+) exp [-Z p ^{r)\ + ^ /3 V 

7'=_ff+,S+ 



sinh (py\V r ip*(r)\) 
Py|V r ^(r)| 



(17) 



Note that /3*(r) + J2y=H+ B+Py( r ) = !• F rom Eqs. [TBI and [JJHwe can interpret that /3*(r) 
and pi(r) denote the probability of finding a charge and a dipole resulting from adsorption 
of counterion of type 7, respectively, at location r. 

Invoking the relation between the Bjerrum length and local dielectric function, 
I b {Pp,ip*} = e 2 /e {p*, ip*} k B T , we obtain 



{p;,r} 



1 + 4vr/ Bo 



PyP y (r) P ;(r)L[p y |V r ^(r) 

y=ij+,s+ 



+P*P S 



1 - 



Ppo 



L\p s \V r r{r) 



(18) 



It is worth noting that the functional form of the field dependent dielectric function in Eq. 



is qua. 



Booth 



31 



itatively similar to the expression derived in the classic papers by Onsager[29| and 
321 ] . Noting that L(x) decreases with an increase in x and sinhx/x (appearing 



in p cf. Eq. [161) increases with an increase in x, the local field can lead to dielectric 
decrement or increment depending on the relative strength of these counteracting effects. 
In other words, the first effect (also arising in the Langevin-Debye model 28j and known as 
the dielectric saturation effect) tends to lower the dielectric function. On the other hand, 
the second effect resulting from the counterion adsorption tend to increase it. 
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An important insight into the effect of electric dipoles on the thermodynamics of the 
polyelectrolyte brushes is obtained if one considers the weak coupling limit for the electric 
dipoles (i.e., the limit of weak dipoles or weak local electric fields) so that py\V r ip* (r)\ — > 0, 
p s |V r '0*( r )| — >* 0- in these limits we can either use the approximation In [sinx/x] — — x 2 /6 
in the expression for the partition function or L(x) — > x/3 for x — > in Eq. [HI Using either 
of the approximations, it can be shown that the electric dipoles renormalize the Bjerrum 
length of the medium by the relation 



1 _ 1 1 4vr 

Ib{p* p ,iJ>*} = h{p*} ' h~ + T 



-y'=H+,B+ 



Ppo 



(19) 



This allows us to write e(r) 



r), where 4> p {r) 



p*(r)/p po and a (r) 



1 — p*(r)/p po are the volume fractions of the monomers and solvent molecules, respectively. 
Also, e p and e s are the dielectric constants of the pure components, given by e p /e Q — 1 = 
47tZbo J2 7 '=h+,b+ (/V#/) Ppo/3 and by e s /e - 1 = 4Trl Bo p 2 s p so /3. In other words, the linear 
mixing rule based on the volume fractions can be used to estimate the local dielectric function 
of the medium in the weak-coupling limit. 

Furthermore, by optimizing the Hamiltonian given by Eq. H] with respect to the varia- 
tional parameter (3 H +, we obtain 



log 



10 



f5 H + Jdrexp [-Z H +ip*(r)] 
1 " Ey N 
1 



n 



H+ 



2.303nN 



drp* p (r) 



/?*(r) 



1 " Ey N 



Ph+( t 



where 7' = H + , B + . In deriving Eq. [201 we have used the fact that nL 



H+ 



n 



total 
H+ 



(20) 



-!3 H +nN. 



Also, we have defined pK^+ = — log 10 K H +. The last term in Eq. [20] captures the effects of 
finite monomer density and inhomogenity on the charge regulation in the brushes. This term 
can be interpretted as the shift in pK H+ and depends on the local environment. However, 
sign of the shift depends on the relative contributions of the charges and dipole. 

A similar equation is obtained for /3b+, which can be written by replacing H + by B + in 
Eq. [201 For the binding of the counterions from the salt, we define pK B + = — log 10 K B +. Eq. 
[20] equates the chemical potentials of the counterions in the "free" state and the "adsorbed" 
state. In the next section, we present the details about the treatment of grafted ends in the 
theory. 
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B. Treatment of grafted ends: uniformly grafted brush 

The functional derivative of the chain partition function with an end grafted at R a (0) 
r a on the substrate allows us to write Eq. [9] as 

J Q N dtq ra (v,t)q(v,N-t) 



P *> [ } ^ fdvq ra (v,t)q(v,N-t) ' 1 } 



where g(r, N — t) satisfies 



dq(r,t f ) 



2 



<)r , QM (22) 

with the condition g(r, 0) = 1 for t' — N — t — 0. Similarly, q Ta (r,t) satisfies the same 
equation but with the initial condition q ra (r, 0) = S(r — r 



Noting that the denominator in Eq. [21] is independent of t, Muller[46j has developed an 
efficient algorithm to solve these equations in three dimensional space. Following Muller, we 
choose t = for the denominator in Eq. [5T] and use the initial conditions for q and q to 
write Eq. [2T1 as 

p* p {r) = / dtq(r,t)q(r,N-t) (23) 
Jo 

where 

^g(r a ,iV) J q{r a ,N) 

Here, J d?2d is the integral over the grafting points. The later transformation to go from a 
discrete sum over all the chains to the quadrature over the grafting points is strictly valid in 
the thermodynamic limit so that n — > oo, v4(=area of the grafting plane) — > oo and a = n/A 
is finite. Due to the linear relation between q and q Ta , q also satisfies Eq. [221 Di rt with the 
initial condition 

«'•■» - rak « 25 » 

where we have written r a = {x a ,y a , 5} in Cartesian co-ordinates so that 5 — » is the 
parameter used to put the anchoring point at an infinitesimal distance above the substrate. 

The above transformation also allows us to write the sum over the chain partition func- 
tions in the Hamiltonian (cf. Eq. H]) in a manageable form, given by 

n „ 

J2^Q pa {w p } ~ a dr 2d \n[q({r 2d ,5},N)} (26) 

a=l J 

= aA\n[q(5,N)] (27) 
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The last expression is valid in the special case when q(r, N) is a function of the distance 
from the substrate i.e., in the absence of lateral inhomogeneities. We call this particular 
case the one-dimensional brush because the densities of different components are dependent 
only on the distance from the substrate. This is the case for high grafting densities in good 
solvents. 



For a polyelectrolyte brush in equilibrium with a solution containing monovalent salt, 
we have to equate the chemical potentials of all the components that can be exchanged 
between the bulk solution and the brush region. Also, at equilibrium, the osmotic pressure 
must be the same everywhere. However, due to the incompressibility constraint, equating 
the chemical potentials of the salt ions and solvent molecules is sufficient to define the 
equilibrium state of the system. Furthermore, it is assumed that the densities of different 
components are known in the bulk solution and we use them as parameters in the study 
here. 

Using the thermodynamic relation between the free energy and chemical potentials in the 
canonical ensemble (i.e., fij = (dF/drij)^), we get 



for the "free" ions, 7 = H + , B + , A~ so that n 7 = for 7 = H + , B + . In writing the second 
equation on the right in Eq. [28], we have used Eq. [13] and the boundary condition that 
the density of the "free" ions is given by p 7 (oo) in the bulk (designated by r — > 00). Also, 
■0*(oo) is the value of ip*(r) in the bulk. The values of /3y are determined by using Eq. 
|2"81 along with Eq. [2D] An important insight into the effect of salt concentration on the 
charge regulation can be obtained by taking the limit of p*(r) = (i.e., the limit of vanishing 
grafting density) in Eq. [20] Taking the limit of zero monomer density in Eq. [20] and using 
Eq. EH we can write 



C. Free energy with respect to the salty solution 



In 



n. 



lnp 7 (oo) + Z 7 '0*(oo) 



(28) 



f drexp [— Z^ip*(r)} 



1 



(29) 



1 + 10P H -P K 



H+ 



e -Z + ^*(oo)/2.303 + p B+ (oo)10P K 



B+ 




(30) 
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where we have defined pH of the bulk as pH = — log 10 [p#+(oo)] and Z+ = \Z H + \ = \Z B + \ is 
the valency of counterions (= 1 for the monovalent ions considered in this work). p B +(oo) 
is the concentration of the cations coming from the added salt in the bulk. 

Eqs. [29] and [30] capture the physics of competitive counterion adsorption on the poly- 
elecrolyte chains in the limit of vanishing monomer density. Eqs. [15] and [16] capture the 
same in an inhomogeneous medium at finite monomer density. In the next section, we'll 
demonstrate that /3*(r), which determines the charge regulation in polyelectrolyte brushes, 
is highly sensitive to the competitive counterion adsorption captured by Eqs. [inland [201 In 
the particular case, when the counterions from the added salt don't adsorb on the polyelec- 
trolyte chains, pK B+ -oo and (3 B+ -> and (3 H + = 1/ 1 + loPH-pK^+e-^^ 00 )/ 2 - 303 . 
For this case, the degree of adsorption of the counterions from the salt (f3 B + = 0) and pro- 
tons (Ph+) become independent of the bulk salt concentration. However, in general, /3h+ 
decreases and (3 B + increases with an increase in the bulk salt concentration. Furthermore, 
in general, the dipoles formed as a result of the adsorption of counterions from the salt have 
higher moments in comparison to those formed by the adsorption of protons. In the next 
section, this difference in the dipole moments is shown to play an important role in affecting 
the dielectric function within a polyelectrolyte brush. 

Similarly, the chemical potential of solvent molecules is given by 



tUssPso 



- In 47r + In 



ft* 



w ss p s 



-ln47r + lnp so + <(oo). (31) 



f drexp [— w*(r)] 

Taking the homogeneous bulk solution as the reference frame, we can write the free energy 
at the saddle-point as F* = F* e j + AF*, where F* e ^ = Y1 1 =h+ b+ a- s AW 111 this work, we 
are interested in the effect of inhomogeneities on the properties of polyelectrolyte brushes 
and AF* is the contribution to the free energy arising from the inhomogeneous distributons 
of fields. Because the free energy is defined with in an arbitrary constant, we subtract a 
constant to ensure that AF* = in the absence of inhomogeneities. Explicitly, AF* is given 



by 

AF F se ij -\- Fcfoemical -\- F w -\- Fj ons -|- F so \ ven t -)- F e -)- Fp iy : 

where F se if/k B T = [w pp p po /2 — ln47r]niV and 

F chem i C ai = nN \(3 B + In K B + - (1 - H +) In K H +] 



(32) 



+nN 



^/3yln/3 y + l-^/3 v )ln 

.7' \ V 



Eft 



(33) 
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where 7' = H + , B + and 



solvent 



poly 



X P s J dr P ;(r)p*(r) 



7=_H"+,_B+,A- 



<(r)ln 



p*(r) In 



PsW 



_ P70 



P 7 W +P 70 



(34) 
(35) 

(36) 



^ Z 7P ;(r)^(r)+^(r) P ;(r) 

7=_H"+,_B+,/l- 

[ sinh ( Ps |V r y;*(r)|) " 
arpJrjm — - — — , . ... — 
L (Ps|V r ^(r)|) . 

a J dr 2d \n [q({r 2d , 6}, N)} - 1 J dr P ;(r)w;(r) 



8tt/ 



Bo 



where p*(r) = p so [l - p*(r)/p p J. 



(37) 
(38) 



D. Numerical methods 

We have solved the non-linear set of equations assuming lateral homogeneity in the planar 
polyelectrolyte brush. For the brush in equilibrium with the solution bath, the above equa- 
tions are written in dimensionless form and the direction prependicular to the grafting plane 
is assigned as the z axis in Cartesian coordinates. All of the quantities having dimensions 
of length are made dimensionless by dividing them by R go = (Nl 2 /6) 1 / 2 : . Numerical calcu- 
lations ensuring that the electric field far from the brush region approaches zero, demand a 
box length as high as 100-R 9O , especially in low salt conditions such as 1 millimolar (mM). 
For the highest bulk salt concentration of 100 mM, a box length of 20R go is found to be 
sufficient. However, the monomer density decays to zero within a distance of 10 — 15R go from 
the grafting surface for all the cases investigated in this work. Numerical results presented 
in this work have been obtained by taking the maximum box length of 100i? 9O with 512 grid 
points. 

The modified diffusion equation represented by Eq. [22] and the Poisson-Boltzmann equa- 
have been solved by using implicit-explicit scheme known as the extrapolated 



tion (Eq. QZ 
gear method 



47J, 



481 ] . In order to use the extrapolated gear method, we rewrite Eq. IT21 in 
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the form 



dip*(r) 
di 



h{p;,r} 



V r ^(r) 



+ 47r P :(r) 



(39) 



where i is a fictitious time and p*(r) = ^2 1= h+ b+ a- ^iP*^) + Z p f3*(r)p*(r) is the local 
charge density. The steady state solution of Eq. EH] is Eq. [TJJ Before implementing the 
extrapolated gear method in one dimension, we rewrite Eq. [3H] as 



di 2 dz 2 dz 



1 



A, 



Ib{p* p ,V} 



dijj*(z) 
dz 



+ ^p* e {z) 



(40) 



where A r 



= 1/[Ib {p»)V'*}]mm is chosen as the maximum value of inverse of the local 
Bjerrum length at each time step. In the extrapolated gear method, the first term on the 



right hand side of Eq. HO] is treated implicitly and the rest is treated explicit 



particular choice of X max is motivated by the studies on the phase field models[48j], which 



y. The 



require the solution of similar equations. An explicit scheme is used for the initialization of 
the extrapolated gear method. 

An equation similar to Eq. HQ]is written while solving Eq. [22]with the choice A max = 1/2. 
Furthermore, Crank-Nicholson scheme is used to initialize the gear method for the modified 
diffusion equation. Dirichlet boundary conditions are used for q and q at the substrate i.e., 
q{z = 0,t) = q(z — 0, t) — for all values of t. The grafted ends are displaced to the 
first grid point in the solution so that 5 = 0.19R go for a box length of 100R go with 512 
grid points. Also, steps of 10~ 4 and 10 -3 are used for the time stepping while solving the 
Poisson-Boltzmann and the modified diffusion equations, respectively. 

In the absence of external electric fields, we have chosen ip*(oo) = 0. Starting from an 
initial guess for the fields (w*, ip*) and charge parameters ((3h+, Pb+), we use the extrapolated 
gear method to compute i/;* at the next time step, which is used as a guess in the iterative 
scheme. On the other hand, the guessed and the computed values for w*, (3h+ and (3b+ are 
mixed using the simple mixing scheme 43|] to develop a new guess for the next iteration. 



Random numbers have been used as initial guesses for the fields w* and ip*. However, Eqs. 
29] and [30] are used as initial guesses for (3h+ and Pb+, respectively. The iterative procedure 
is continued until the free energy of the brush (with respect to the solution bath) doesn't 
change within 10~ 8 . 
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IV. RESULTS 

Although the theory is quite general, in the following we have applied it to study a brush 
made of chains bearing acidic groups having pK H + = 4.66, which mimics poly(methacrylic 
acid) The motivation behind this choice is the availability of data for the stability 



constants 



50| of monovalent salts for this system and on-going experiments. For monovalent 



salts such as sodium chloride (NaCl) and lithium chloride (LiCl), it is found 50|, [51J that 
pK Na+ ~ pK Li + = 0.28. Furthermore, we present the results for p H + = 0.035nm = 1.7 
Debye and p s = 0.1003nm = 4.8 Debye mimicking weakly acidic groups such as methacrylic 



acid 44J and water, respectively, at the room temperature. The temperature is fixed in 
the calculations by the choice of Ib — 56nm, which corresponds to room temperature. 
The choice of p s = 0.1003 nm is motivated by the fact that the dielectric constant of 
water is 80 at room temperature. Note that the value of p s chosen for the numerical work 
here is little higher than the actual dipole moment of water (= 1.85 Debye in the gas 
phase )• The origin of this lies in the neglect of hydrogen bonding and induced dipole 
moment effects in the theory. Overcoming this limitation of the theory is an interesting 
direction for future research. Also, we have taken I = 0.3103 nm, which corresponds to 
Pso — Ppo = l// 3 = lg/cm 3 i.e., the density of water at room temperature. Note that the 



choice of / = 0.3103 nm is also close to the Kuhn segment length 



53| of poly(methacrylic 



acid). We have solved the equations at the saddle-point in three dimensional space for 
N = 100 and our preliminary studies show that a grafting density of uRg = 1 for Xpsl 3 = 0.6 
is high enough to ensure that there are no lateral inhomogeneities. Note that the choice of 
Xpsl 3 = 0.6 represents a slightly poor-solvent quality for the polyelectrolyte chain backbone, 
which is definitely the case for poly (methacrylic acid) in an aqueous medium. All the results 
presented below are for iV = 100 and Xpsl 3 — 0.6. 

Using these set of parameters relevant for poly(methacrylic acid) in an aqueous medium, 
we have varied the salt concentration and the pH of the bulk solution. To study the effect 
of adsorption of counterions coming from the salt, we have considered two cases. The first 
corresponds to pK B+ = 0.28 and the second corresponds to pK B+ — > —00, which represents 
the scenario in which the counterions from the salt do not adsorb on the chains. For the first 
case, where counterions from the salt can also adsorb on the acidic chains, we have taken 
Pb+ = 0.125 nm (corresponding to dipole moment of 6 Debye for the sodium carboxylate 
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-COONa) group |52|). 



A. Effects of the bulk salt concentration 



1. Monomer density profiles, charge regulation and the dielectric function 

In Fig. |2] we present the results obtained for the polyacidic brushes when the bulk 
salt concentration,^ = p#+(oo), is varied. These results correspond to bulk pH = 5 and 
pK B+ — > —oo so that counterions coming from the salt don't adsorb on the chains. 

Monomer density profiles (Fig. [2(a)) show a depletion zone near the substrate, which is 
an outcome of using Dirichlet boundary conditions for the chain propagator at the substrate 
and a delta function as the initial condition. These boundary conditions correspond to a 
non-adsorbing and highly repulsive substrate. Also, the monomer density profiles approach 
zero around z/R go = 6 — 10, depending on the bulk salt concentration. Furthermore, the 
maxima in the monomer density profiles decreases in magnitude with an increase in the bulk 
salt concentration. In other words, the height of the brush as manifest in the total extent 
of the chains increases with an increase in the bulk salt concentration. This is an outcome 
of an increase in the local charge inside the brush as the bulk salt concentration increases, 
which is reflected in Fig. |2jc). 

The monomer density profiles affect the dielectric function in a significant manner as 
shown in Fig. [2jb). First, it is worth noting that the dielectric function closely follows 
the monomer density and a dielectric decrement is observed inside the brush region. The 
dielectric decrement is an outcome of the fact that the acidic groups on the chains are less 
polar than the solvent (ph+ < Ps)- The actual magnitude of the local dielectric function 
depends on the local density, electric field, temperature and the dipole moment of ion-pairs. 
In this work, we have fixed the temperature at room temperature by choosing Ib = 56 nm. 
In the region far from the grafting plane where the monomer density is zero, a dielectric 
constant of 80 is obtained, representing solvent rich areas. Inside the brush, the dielectric 
decrement is as high as 30, and it is seen that the dielectric decrement decreases with an 
increase in the bulk salt concentration. The decrease in the dielectric decrement is a result 
of a decrease in the local monomer density due to an increase in the local charge, as shown 
in Fig. 12(c). 
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FIG. 2: Effect of the bulk salt concentrations on properties of the polyacidic brush in the case 
where counterions from the salt don't adsorb on the chains. These results are obtained for bulk 
pH = 5 and pK B + — > — oo. Other parameters are presented in the main text. Figs. (a),(b),(c) 
and (d) correspond to the monomer density (p*(z)), local dielectric function (e*(z)), probability of 
finding a monomer in a charged state (/3*(z)) and the electrostatic potential (ip*(z)), respectively. 
The saddle-point equations are solved for z max /R go = 100 and the data is cut at z/R go = 15 and 
50 to highlight the important features for the Figs, (a)-(b) and (c)-(d), respectively. 



From Fig. |5] (c) it is found that the probability of finding a monomer in the charged state 
(f3*(z)) is higher in the bulk solution and decreases inside the brush region followed by an 
increase in the depletion zone near the substrate. This is in agreement with the fact that the 
degree of ionization of polyelectroiytes decreases with an increase in monomer density 411] . 
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Also, comparing the results from the calculations for different bulk salt concentrations, it is 
clear that the magnitude of /3*(z) increases with an increase in the bulk salt concentration. 
In other words, extent of charging of the polyelectrolyte chains increases with an increase 
in the bulk salt concentrations. This is a result of the decrease in the electrostatic potential 
with an increase in the bulk salt concentration (Fig. [2] (d)), which leads to a shift in pK H+ 
(Eq. 12U1) . Similar effects of the salt concentration on the charge regulation are observed in 



recent theoretical calculations 



161 ] and experiments 40J . 



The non-monotonic behavior of the probability of finding a monomer in the charged state 
is closely related to the non-monotonic behavior of the electrostatic potential as shown in 
Fig. [2] (d). It is observed that the electrostatic potential is negative everywhere for the 
polyacidic brush and approaches zero as far as zj R go = 50 for the cases where the monomer 
density approaches zero at z/R go ~ 6 — 10 (cf. Fig. [2] (a)). This long-range behavior of 
the electrostatic potential complicates the numerical investigation of the strongly charged 
polyelectrolyte brushes in the absence of salt. We haven't extended the study to salt-free 
strongly charged polyelectrolyte brushes precisely because of this reason. However, we expect 
the results to be qualitatively similar. 

Inside the brush region (i.e., z/R go < 6 — 10 depending on the bulk salt concentration), 
the electrostatic potential closely follows the monomer density profiles (compare Figs. [2] 
(a) and 121(d)). The minima in the electrostatic potential corresponds to the maxima in the 
monomer density profiles and inside the depletion zone, the magnitude of the electrostatic 
potential decreases due to a decrease in the monomer number density. An increase in the 
bulk salt concentration leads to an increase in the local degree of charging ((3*(z)) and the 
electrostatic potential also increases in magnitude, as seen in Fig. 121(d). 



2. Counterion and coion distribution 

Figs. |3](a) and [3](b) show the density profiles for the "free" counterions (i.e., positively 
charged H + and B + ions) and coions (i.e., negatively charged A~ ions), respectively, for 
different concentrations of the salt as set by the bulk solution. Physically, electrostatic 
attraction between the charged monomers and the counterions are responsible for an increase 
in the density of the counterions inside the brush region in comparison with the bulk. 
Similarly, electrostatic repulsion between the charged monomers and co-ions are the reason 
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FIG. 3: Distributions of small ions at pH = 5 and different salt concentrations in the case where 
counterions from the salt don't bind on to the chains. Figs, (a) and (b) show the total counterion 
and co-ion density profiles, respectively. The parameters used for producing these figures are the 
same as in Fig. [2] and the data is cut at z/R go = 50. 



behind a lower number density of co-ions inside the brush region compared to the bulk 
solution as shown in Fig. [3(b). 

The magnitudes of the electrostatic attraction and repulsion for the monomer-counterion 
and monomer-co-ion pairs, respectively, are related to the electrostatic potential profiles 
shown in Fig. 12(d). The non-monotonic behavior of the electrostatic potential manifests in 
the counterion and co-ion density profiles. Furthermore, the strength of these interactions 
is modulated by a change in the bulk salt concentration. For example, an increase in the 
bulk salt concentration leads to an increase in the degree of charging of monomers (/3*(z)) 
and electrostatic potential (ip*(z)), and increases the effect of electrostatic attractions and 
repulsions, as seen in Fig. [3J 

At this point, we comment on an important issue not considered in this work. We note 
that non-trivial effects such as the partitioning of small ions based on the ion solvation en- 
ergy are not taken into account. However, we can infer the effects of solvation energy by 
the fact that solvation energy of small ions (within Born's theory 44J] ) is inversely propor- 
tional to the dielectric constant. This, in turn, means that the ions prefer regions of higher 
dielectric functions over regions of low dielectric functions. In the case of the co-ions in- 
side a polyelectrolyte brush showing dielectric decrement, this effect works in synergy with 
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the electrostatic repulsion between the monomers and the ions. In contrast, in the case 
of the counterions, the solvation effect counteracts the electrostatic attraction between the 
monomers and the coions. Depending on the relative strength of these counteracting effects, 
the counterion density may show a non-monotonic behavior when the solvation effects are 
taken into account. However, treatment of the solvation energy of charges on the chains 
complicates the analysis and this important direction of research is reserved for future work. 

B. Effects of the competitive counterion adsorption 

Fig. H] shows different characteristics of the poly(acid) brushes for pK B+ = 0.28 so 
that counterions from the added salt can also adsorb on the polyelectrolyte chains. From 
Figs. [2] and HI it is found that the monomer density profiles and the dielectric function 
profiles are qualitatively similar in the two cases of counterion adsorption. However, there 
are quantitative differences and in particular, the charge regulation is affected in a significant 
manner as seen in Fig. 11(c) . For example, with an increase in the bulk salt concentration, 
/3*(z) in the bulk solution, first increases (for c s < 10 mM) and then decreases (for c s = 
50, 100 mM). However, inside the brush region, fl*(z) increases with an increase in the salt 
concentration, as seen in Fig. [2](c) also. The increase in the fl*(z) with an increase in the salt 
concentration should be interpretted as a screening effect, which modifies the electrostatic 
potential (cf. Fig. 11(d) ) and in turn, the extent of counterion adsorption. This effect 
depends on the concentration of "free" ions (see Figs. [5(a) and [5(b)). However, in the case 
where counterions from the salt can also adsorb, Eqs. [2H]and [313 which are strictly valid at 
zero monomer density, reveal that an increase in the salt concentration at a fixed pH leads 
to the adsorption of more counterions from the salt due to enhanced dissociation of protons. 
This effect is observed in Fig. 3(c) for finite monomer density at c s = 50, 100 mM, where 
the bulk value of 0*(z) decreases with an increase in the bulk salt concentration. 

In other words, an increase in the bulk salt concentrations causes the counterion adsorp- 
tion equilibrium to shift toward an increase in the adsorption of the counterions from the 
salt at the expense of dissociation of protons. Due to the fact that ion-pairs formed due 
to the adsorption of counterions from the salt are more polar than the native acid groups, 
the dielectric function is also affected by the shift in the counterion adsorption equilibrium. 
This can be observed by comparing Figs. [2(b) and 11(b). For example, at c s = 100 mM the 
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FIG. 4: Effect of the bulk salt concentrations on properties of the polyacidic brush in the case 
where counterions from the salt can also bind for the bulk pH = 5 and pKe+ = 0.28. All the other 
parameters are the same as in Fig. [2j Figs. (a),(b),(c) and (d) correspond to the monomer density 
(Pp(z)), local dielectric function (e*(z)), probability of finding a monomer in a charged state (f3*(z)) 
and the electrostatic potential (ip*(z)), respectively. Like Fig. [2j the data is cut at z/R go = 15 
and 50 to highlight the important features for the Figs, (a)-(b) and (c)-(d), respectively. 

lowest local dielectric function inside the brush is 65 for pK B + = 0.28 in comparison to the 
value of ~ 50 for pK B+ — > — oo. In the next section, we demonstrate that the shift in the 
counterion adsorption equilibrium from protons toward salt ions can also lead to dielectric 
increment just by modulating pH of the bulk solution. 
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FIG. 5: Distributions of the counterions (a) and co-ions (b) for pH = 5,pK B + = 0.28 and different 
salt concentrations in the bulk. Monomer densities and other properties of these brushes are shown 
in Fig. |U 



It is well-known that an increase in pH of the bulk solution leads to an increase in 
the charge on poly(acid) chains. This, in turn, means that there are more sites for the 
adsorption of counterions from the salt. In Fig. El we present the results of calculations 
where dielectric increment is observed due to the replacement of less polar acidic groups by 
more polar ion-pairs. 

From Figs. EJa) and E^c), it is clear that both the height of poly (acid) brushes and local 
charge increase with an increase in pH of the bulk solution. Furthermore, the electrostatic 
potential increases in magnitude as the bulk pH increases. However, the dielectric function 
inside the brush region is greater than that of the bulk for pH = 6, 7 in contrast to being 
lower for the other pH values. As mentioned earlier, this is a direct outcome of salt ion 
adsorption onto the chains. 

An increase in the local charge with an increase in the bulk pH affects the counterion 
and co-ion density profiles, as shown in Fig. [7J An increase in the local charge leads to 
an increase in the electrostatic potential, which manifests as an increase in the counterion 
density profiles and as a decrease in the co-ion density profiles. 



C. Effects of the pH of the bulk solution 
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FIG. 6: Effect of the bulk pH on the properties of poly(acid) brushes at a bulk salt concentration 
= 50 mM and pK B + = 0.28. All the other parameters are the same as in Fig. [2j Figs. (a),(b),(c) 
and (d) correspond to the monomer density (pp(z)), local dielectric function (e*(z)), probability of 
finding a monomer in a charged state (f3*(z)) and the electrostatic potential (ip*(z)), respectively. 
Like Fig. [2j the data is cut at z/R go = 15 and 50 to highlight the important features for the Figs. 
(a)-(b) and (c)-(d), respectively. 

D. Importance of non-linear effects 

To highlight the importance of non-linear effects such as the dependence of the dielectric 
function on the local field and charge, we have compared our results with two approaches 
based on the assumption that the weak-coupling limit is valid independent of the coupling 
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FIG. 7: Ion distributions for different pH values in the bulk at a salt concentration = 50 mM in 
the case where counterions from the salt also bind (pK B + = 0.28). Figs, (a) and (b) represents 
the counterion and co-ion density profiles, respectively. 

strength, determined by p 7 /|V r ^*(r)|. In particular, we have compared (Fig. E]) the di- 
electric function predicted using the non-linear theory described herein with the behaviors 
expected from linear mixing rule for the dielectric function, which is shown to be valid in 
the weak coupling limit (cf. Eq. H~9]) . Also, to highlight the importance of non-linear effects 
on the charge regulation and, in turn, on the dielectric function, we have compared the 
results obtained by calculating (3 H + and using either Eq. [2D] (from the full non-linear 
theory) or Eqs. [29] - [30] strictly valid in the limit of zero monomer density. These two 
approaches for the computation of dielectric function based on different estimates for (3h+ 
and Pb+ m the weak coupling limit are termed as WCLC (weak coupling limit with charge 
regulation) and WCLD (weak coupling limit in dilute regime), respectively. Note that the 
computations of the dielectric function in the WCLC and WCLD methods are done after 
solving the SCFT equations for the full non-linear theory. In case, the non-linear effects are 
insignificant, the results obtained from the three methods (i.e., WCLC, WCLD and the full 
non-linear theory) should be identical. 

As seen from Fig. [8] it is clear that the non-linear effects are not important when the 
degree of dissociation is low such as when pH= 3 (cf. Figs. |S{a) and Etc)), as expected. In 
this case, non-linear effects on the charge regulation and, in turn, on the dielectric function 
are insignificant. In contrast, at a moderate degree of dissociation (e.g., for pH = 5 in 
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FIG. 8: Significance of the non-linear effects on the dielectric function is demonstrated here. WCLC 
and WCLD represent two different ways of computing the dielectric function in the weak coupling 
limit. For the details, see the main text. The other plots for different pH values are the same as 
in Fig. Mjo). 

Fig. lEfb)), there are noticeable differences between the dielectric functions predicted by 
the different approaches. In particular, comparing the dielectric functions obtained from 
the WCLC and WCLD methods in Fig. Efb), it is evident that the non-linear effects on 
charge regulation (arising from finite monomer density, cf. Eq. |20|) play an important role. 
Furthermore, additional contributions arising from non-linear dependence of the dielectric 
function on the local field also play a significant role. These effects lead to a maximum 
difference of ~ 5 in the dielectric function, as seen in Fig. [8(b). The non- linear effects are 
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the most significant in cases where the charged groups are fully dissociated such as when 
the brush is in a solution at pH = 7, which is shown in Fig. E(c). In this case, the non- 
linear theory predicts a dielectric increment inside the brush, which is in striking contrast 
to the dielectric decrement predicted by the WCLC and WCLD methods. The origin of 
this discrepancy lies in the neglect of the dependence of local dipolar density on 
the local field in the WCLC and WCLD methods. The dipolar density increases with an 
increase in the coupling strength (p 7 /|V r , i/'*(r)|) as seen from Eq. [TBJ However, the non- 
linear effects are not significant far from the brush regime where pylVrV^r)! — >■ and the 
results obtained using the three different methods are the same, as expected. 

V. CONCLUSIONS 

We have studied planar polyelectrolyte brushes (made of end-tethered poly(acid) chains) 
in equilibrium with an electrolyte solution using field theory. We have focused on the 
quantitative description of the charge regulation and local dielectric function. Although the 
theory is quite general, in this work we have studied the effects of salt concentration and pH 
of the bulk assuming a laterally homogeneous brush. In addition, the effects of competitive 
counterion adsorption are studied by allowing the salt ions to form ion-pairs by adsorption 
onto the charged monomers. 

The dipole moment of the ion-pairs is shown to significantly affect the dielectric function. 
For the poly(acid) chains bearing groups less polar than the solvent, the dielectric function 
inside the brush region is predicted to be lower than the bulk solution in the absence of 
adsorption of any salt ions. However, the formation of ion-pairs, generally, more polar 
than the solvent molecules, is shown to increase the dielectric function inside the brush in 
comparison to the bulk solution. Comparison of the theory taking into account dependence 
of the dielectric function on the local field and charge regulation (non-linear effects) with 
other approaches treating dielectric function by the linear mixing rules reveals that the non- 
linear effects are significant and must be taken into account in order to predict qualitatively 
correct behavior (such as the dielectric increment inside the brush) in charged systems. 
Furthermore, it is shown that local charge inside the brush region shows non-monotonic 
behavior. The local charge inside the brush region increases with an increase in the bulk 
salt concentration. However, the counterion adsorption equilibrium shifts from protons to 
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the salt ions with an increase in the salt concentration. Counterion and co-ion densities are 
predicted to be higher and lower, respectively, in comparison with their bulk value. 

We also comment on an important assumption used in the theory. We have ignored 
the effects of "induced" dipole moments, which have been shown to enhance net dipole 



moment 



291 ] . Ignoring this induction effect is a reasonable first step before developing a 



more comprehensive theory to treat the dielectric response of flexible macromolecules. In 
the future, we plan to overcome this limitation of the theory by extending the field theoretical 
method developed here. 

Finall y i t is understood that the dielectric function has an intricate relation with the 



solvation 
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U457f of charges. Classic work by Born[44| directly relates the solubility of 
charged molecules into solvents of different dielectric constants and predicts the partitioning 
of salt ions based on their valency and radii. In the literature it has been shown that the 



solvation of ions can lead to non-trivial and counter-intutive results. For examp 
induced stabilization of the ordered morphologies in the mixtures of small molecules 



e, salt 
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is shown to be an outcome of the asymmetric 
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63] solvation of the cations and anions. 



The theoretical treatment presented here employs the saddle-point approximation for point 
ions and dipoles. Thus, the theory in the current form is unable to capture effects such as 
asymmetric solvation of ions based on their radii. However, the theory can be extended to 



finite size ions by using a soft primitive model for charges 
issue in a future publication. 



6J-66|. We plan to address this 
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APPENDIX A : Counterion adsorption, dipolar interactions and the partition 

function for polyelectrolyte brushes 

Here we present the partition function for the polyelectrolyte brushes in the presence of 
polar solvent molecules and small monovalent ions. Special attention is paid to the counte- 
rion adsorption and dipolar interactions. We start from a molecular description taking into 
account the dipolar interactions along with other short range interactions (dispersion and 
van der Waals). In order to cast theory in a field theoretical language, we rewrite the parti- 
tion function in a form suitable for mathematical transformations. In this work, we consider 
polyelectrolyte brushes made of n mono-disperse end-tethered chains, each containing N 
Kuhn segments. 



Following Edward's work[42j, we represent a polyelectrolyte chain in the brush by a 
continuous curve R(£) of length Nl, I being the Kuhn segment length and t is the arc 
variable representing any point on the curve lying in the range (0, N). To keep track of 
different chains, subscript a is used so that t a represents the contour variable along the 
backbone of a th chain. Also, the position vector for a particular segment ,t a , is written as 

In the following, subscripts p, s and 7 are used to represent monomers, solvent molecules 
and the small ions, respectively. Three different kinds of small ions are considered here and 
unless specified, 7 = H + , B + and A~ represents protons, cations from the salt and anions, 
respectively, as described in the main text. In this work, we study negatively charged (or 
polyacidic chains) chains and specificity of the cations is taken into account to study the 
effects of different binding energies of the cations. For the treatment shown below, we assume 
that local incompressiblity condition is satisfied. Treating all the small ions as point-like, 
the incompressiblity condition allows us to write the total volume as Q = nN/ p po + n s / p so 
where p po and p so are the bulk monomer and solvent densities, respectively. Also, we use the 
notation / 3 = l/p po and represents the position vector of k th small molecule like solvent 
molecules and small ions. 

Furthermore, counterion adsorption on the polyelectrolyte chains is taken into account 
using a two-state model. Segments along the chains can be either in charged or in uncharged 
state. To describe the two states, another arc length variable, a (t a ), is introduced, which 
enumerates the state of charging of the segment, t a on a th chain. For the analysis here, 
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Qa(t a ) = means t a is a neutral site and 9 a (t a ) = 1 represents a fully charged site along the 
backbone. Like the average over all of the possible conformations in the theories of neutral 
polymers, we need to average over all of the possible charge distributions along the chains. 
We represent the average over 9 a (t a ) by symbol J2{e a } ((' ' '))> which expicitly means 

]Jdt a D [9 a {t a )\ (■ • -)P [9 a (t a )] T [9 a (t a )} . (A-l) 

Here, P [9 a ] is the probability distribution function for the variable 9 a , which must satisfy 
the relation J D [9 a (t a )} P [9 a (t a )] = 1. Also, T [9 a (t a )] is the number of indistinguishable 
ways in which 9 a can be distributed among nN sites for a fixed number of charged sites. 
T [0a (£<*)] takes into account the entropy of distribution of charged sites. 

Like the segments, the counterions are also divided into two sets. One set of counterions 
is "free" to explore the whole space and has translational degrees of freedom. The other 
set is "adsorbed" on the backbone (ion-pairs) and behave as electric dipoles. Number of 
counterions in "free" and "adsorbed" states are taken to be and n", respectively, for 
7 = H + , B + . In this work, we treat the ion-pairs as point electric dipoles for the development 
of an understanding of dielectric function. In the following, the dipole moment along the a th 
chain backbone is written as by vector, p a (t a ). Similarly, represents the dipole moment 
of the k th solvent molecule. 

Using the notations described above, the partition function (Z) for a polyelectrolyte brush 
can be written as 

n I „ n N „ n 1 i+n s +n A - „ s 

n D[R a] e ( / n n / n <]„ , . n / 

a=i {e a (t a )} V a=it a =o J [[ 7 'n 7 An A -\n s \ j=1 J k=l 

exp [ Hq -[Rq,]- H w {Rq,, R a '} H C p {R Q , p a , r^} Hpp {R Q , 

Pa, R a' ; Pa'} 



-H cc {r v r' J }-E{9 a }]\[5 



Ppo 



(A-2) 



where 7' = H + , B + . 

The Hamiltonian in Eq. IA-21 is written by taking into account the contributions coming 
from the chain connectivity (given by H in Eq. IA-31 below), the short ranged dispersion 
interactions (represented by H w in Eq. IA-41) and the long range electrostatic interac- 
tions between the charged species (written as H pp , H cp and H cc above, which correspond 
to the dipole-dipole, charge-dipole and charge-charge interactions, respectively). For conve- 
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nience in writing, in the following we have suppressed the explicit functional dependence of 
Hq, H w , H pp , H cp and H cc . 

Explicitly, contributions from the chain connectivity are given by: 

a=l ' 







which is the well-known Wiener measure for a flexible polymer chain. Furthermore, 
H w {R Q , R a '} takes into account the energetic contributions coming from short range dis- 
persion interactions between segments of chains indexed as a and a' located at R Q and 
Rq/, respectively. Following Edwards's formulation 42J for a flexible chain, we model these 
interactions by delta functional/point interactions, which allows us to write 

H w = ~ J dr[wpppl(r)+w ss pl(r) + 2w ps p s (r)p s (r)] (A-4) 

Here, w pp , w ss and w ps are the well-known excluded volume parameters describing the 
strength of interactions between p — p, s — s and p — s pairs, respectively. Also, p p (r) and 
p s (r) represent the microscopic number density of the monomers and the solvent molecules, 
respectively, at a certain location r defined as 

r-jV 



p p {r) = y~] I dt a S [r - R a (t a )} 



(A-5) 



p.(r) = ^5[r-r fc ] (A-6) 

k=l 

Electrostatic contributions to the Hamiltonian arising from the dipole-dipole interactions 
can be written by considering finite/real dipoles first and then systematically approaching 
the limit of point dipoles so that the dipole length goes to zero and the charge increases 
keeping the dipole moment finite. For real dipoles, we can explicitly add the Coulomb 
interaction energies between the charges at the ends of dipoles on the a th polyelectrolyte 
chain and solvent molecules. The magnitude of the charges (in units of electronic charge, 
e) at the end of dipoles on the chains and solvent molecules are written as Z p k and Z s k, 
respectively so that k — ± and Z p+ = —Z P ^,Z S+ = —Z s _. For the dipoles located on 
the polyelectrolyte chains, we consider their centers to lie on the chain backbone, and the 
dipolar axes along the vector p a {ta)- Summing over all the dipoles, we obtain the dipole- 
dipole interaction energy as 
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H, 



VP 



Ibo 
~2~ 

+ 



dr / dr' dp dp' 



+ 



7 2 



+ 



7 2 

Z P- 



r - r' + 0.5p - 0.5p'| r - r' - 0.5p + 0.5p'| 



Zp+Zp- 



r — r' + 0.5p + 0.5p'| |r — r' — 0.5p — 0.5p 

zL z? 



7|| Pp(r,p)p P (r / ,p / 



+ 



|r - r' + 0.5p - 0.5p'| |r - r ; - 0.5p + 0.5p' 

Z s +Z s - 



Z s +Z s - 



|r - r ; + 0.5p + 0.5p'| |r - r' - 0.5p - 0.5p 



7|| Ps(r,p)Ps(r',p') 



+2 



Zp+Z s ^ 



Z p -Z s - 



|r - r' + 0.5p - 0.5p'| |r - r' - 0.5p + 0.5p'| 

T\ \ Pp( r 'P)ps( r 'iP' 



Z p +Z s _ 



Zp-Z sA 



(A-7) 



|r - r' + 0.5p + 0.5p'| |r - r' - 0.5p - 0.5p' 

where Ib — e<2 I^-o^b^ is the Bjerrum length in vacuum and e Q is the permittivity of the 
vacuum. Also, ksT is the Boltzmann constant times the temperature. Also, in Eq. \A-7\ we 
have defined the microscopic number density of dipoles on the chains and solvent molecules 
by p p (r, p) and p s (r, p), respectively. Physically, these functions describe the number of 
dipoles with their centers at a certain location r with dipole moments given by p. Formally, 
these are defined as 

n „ N 

p p (r, p) = V / dt a 6[r- R a {t a )\ 8[p- p a {t a )\ [1 - 9 a {t a )\ (A-8) 
«=i Jo 

p s (r, p) = 5 t r - r d 6 [P - P*l ( A " 9 ) 

k=l 

Repeating the same arguments as above, we can write the interaction energy between the 
charges and the dipoles as 



H, 



<v 



I Bo I dr I dr' / dp' 



^Z 7 p 7 (r) + Z p p pe (r) 



Z P +p P (r',p' 



r — r' — 0.5p'| 



+ 



Zp-p P (r', p' 



+ 



Z s+ p s (r',p') Z s _p s (r',p' 



(A-10) 



|r-r' + 0.5p'| |r-r'-0.5p'| |r - r' + 0.5p'| 
where p 7 (r) represents the local microscopic densities for the point-like ions of type 7 at r, 



defined as 



P7W 



r — r. 



for 7 = iJ + ;J B + ,A-. 



(A-ll) 
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Furthermore, p pe (r) is the contribution to the charge density coming from the polyelectrolyte 
chains, given by 

n r N 



/ dt a 5[r-R a (t a )]6 a (t a ) 

a=l Jo 



(A-12) 



Electrostatic interaction energy among the small ions can be written using Coulomb's 
law as 



^ / r/r / <lr'- 



E 7 ^ 7 p 7 (r) + Z p p pe {r) 


^ 7 Z 7 p 7 (r') + Z p p pe (j') 




r — r' 





(A-13) 



The electrostatic terms, written for "finite" electric dipoles, can be written in a simplified 
form, if we consider the limit of "point" dipoles so that the length of the dipoles approaches 



zero. Also, as mentioned earlier, for the electric dipoles, Z p 



Z p and Z s 



—Z s - = Z s . Approaching the limits of infinitesimal dipole lengths in such a way that the 
dipole moments approach their finite values, we can use multipole expansion and write 



H e = H pp + H cv + H, 



cp 



'fUr 



dr' 



Vr'.-Pa 



r — r 



(A-14) 



where p e (r) = Z 7 p 7 (r) + Z p p pe {v) is the local charge density Also, P ave (r') = 
J dpP(r, p)p, is a vectorial quantity so that P(r, p) = p p (r, p) + p s (r, p) is local dipole den- 
sity at r with the dipole moment p. Such an expression for the electrostatic contributions 
resulting from polarization was proposed by Marcus |67| and Felderhof[68j. Furthermore, 



71] using this 



statistical mechanical studies have been carried out for small molecules [69 
expression for the electrostatics in the free energy. 

Note that the electrostatic terms written above depends on the arc length variable 9 a (t a ). 
This variable also determines the energetic contributions of counterion adsorption on the 
backbone, written as E {6 a } in Eq. IA-21 Noting that the dissociable groups on the chains 
have to dissociate first for the salt ions to adsorb, it is written as 

E {9 a } = {nN - n a H+ ) [p° coo . + p° H+ - p° COOH ] + n a B+ [p° coo - B+ - p° coo - - /4+](A-15) 

where n 7 is the number of counterions of type 7' = H + , B + adsorbed on the backbone. p° is 
the chemical potential (in units of fc^T) for species of type j in infinitely dilute conditions. 
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The differences in the chemical potentials are related to the equilibrium constants of the 



corresponding reactions by the relations 



45| 



Vcoo- + V°h + ~ »°cooh = 2.303pK^ + = -2.303 log 10 K H+ (A-16) 
Vcoo- + V°b+ ~ V°coo-B + = 2.303pK B+ = -2.303 log 10 K B+ (A-17) 

Here, we have defined K H + and K B + as the equilibrium (dissociation) constants for the 
reactions 

- COOH ^ -COO- + H + (A-18) 

-COO-B+ ^ -COO- + B + (A-19) 

respectively. Such a model of counterion adsorption was developed in the classic paper by 
Harris and Rice [72]. 

APPENDIX B : Field theory for polyelectrolyte brushes : annealed charge 

distribution 

The probability distribution, P, needs to be determined self-consistently by the mini- 
mization of the free energy and must satisfy the relation J D [0 a (t a )] P = !• m this 
work, we take a variational ansatz for P and write it as 



\yY=H+,B+ J \ -y'=H+,B+ J 

so that n't, = (3ynN for 7' = B + . Mathematically, (3h+ and /?b+ are the variational pa- 
rameters, which will be determined by minimization of the free energy. Physically, and 
Pb+ correspond to the fraction of sites on the chains occupied by H + and B + , respectively. 
Furthermore, using Eq. IB- II for the charge distribution 

nN\ 

T[ ^ (ta)] = n- H+ \n a B+ \(nN-n a H+ -n a B+ )\ (B " 2) 



Such a distribution is called "annealed" distribution in the literature [73l475| . The use of Eq. 
IB- II for the probability distribution the charged sites on the chain backones simplifies the 
subsequent analysis. Using Eq. IB-14 we can write Eq. I A- 2 1 in field theoretic form. We start 
from the electrostatic contribution to the partition function, which is given by Eq. IA-141 
and use the Hubbard-Statonovich transformation 
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exp [-H e 



where 



D 



exp 



dr |p e (r) - V r .Pave(r)}^(r) + 



Bo 



dr^(r)V^/>(r)) 



(B-3) 



/ D[ % 






b] exp 


[ 1 / 







dr^(r)V^(r)) 



(B-4) 



Using this transformation and Eq. IB- 1 ^ we can integrate over the orientations of the dipoles 
analytically and evaluate the average over 9 a . Writing the magnitude of the dipole moments 
of solvent molecules as p s and ion-pairs as py for 7' = H + , B + , the partition function given 
by Eq. IA-21 becomes 

/n n 7 n s if 

n ^[ r j / n n ^ n ^ / ° m ^ t-^ ^ - f « { r -> r ^ 
a=l J 7 j=l k=l " J 

-i J rfr^Z 7 p 7 (r)^(r) -z ^ dr/5p(r)^ p (r) + ^ dr^(r)V^(r)) 

7 
sin (p s |V r ^(r)|) 



+ / drp s (r)ln 



p s |V r ^(r)| 



n* 



(B-5) 



where 



k B T 



n a B + lnK B + - (nN - n a H+ ) In K#+ - In 



n a H+ \n a B+ \{nN — n a H+ — n B+ )._ i 



+ ln 



7i^ + !n£+!n a ! 



- (niV + n a ) ln47r 



(B-6) 



and we have absorbed numerical prefactors of 47r coming from orientational degrees of free- 
dom of the dipoles in the definitions of K#+ and K B +. Also, 

sin (p 7 |V r ^(r)|) 



exp [-# p (r)] = (1 - /3 H + - p B +) exp [-iZ p ip(r)) + ^ /3 7 



(B-7) 



p 7 |V r ^(r)| 

Now, we introduce collective variables corresponding to p p and p s by using the identity 
Z{p P ,Ps} = f D[ Pp ] J D[p s ]Z{p p ,p s }]j6[p p (v)-p p (r)}[p s (r)-p s {r)) (B-8) 

Also, the local incompressibility constraint allows us to rewrite H w given by Eq. IA-4l in terms 



of collective variables as 



[wpppp nN + w ss p so n s } + x P s / drpp(r)p s (r), 



(B-9) 



where we have used f dr Pp (r) = nN, f dr Ps (r) = n s . Also, we have defined a parameter Xps 
by 



Xps tu. 



WppPpo ^sspso 



ps 



(B-10) 



which has the dimensions of volume. 

Now, writing the local constraints (represented by delta functions) in terms of functional 
integrals by the identities 



Y[S[ Pj (r) - Pj (r)] = J D[ Wj ]exp -i J dr Wj (r) { Pj (r) - Pj (r)} 



for j = p,s 

(B-ll) 



and computing trivial functional integrals over Ps (r), we arrive at Eq. [TJ Furthermore, note 
that in the absence of dipolar species (i.e., pj = for j = H + , B + , s in the above treatment), 
•consistent field theory for "annealed" charge distributions documented 



we recover the sel: 
in literature 
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